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The relativistic recoil contributions to the Uehling corrections are revisited. We consider a con- 
troversy in recent calculations, which are based on different approaches including Breit-type and 
Grotch-type calculations. 

We have found that calculations of those authors were in fact done in different gauges and in some 
of those gauges contributions to retardation and two-photon-exchange effects were missed. We have 
evaluated such effects and obtained a consistent result. 

We present a correct expression for the Grotch-type approach which produces a correct gauge- 
invariant result. 

We also consider a finite-nuclear-size correction for the Uehhng term. 

The results are presented for muonic hydrogen and deuterium atoms and for muonic helium-3 
and helium-4 ions. 
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I. INTRODUCTION 

A recent experiment performed at PSI on muonic- 
hydrogen Lamb shift [l| has reported a high-precision re- 
sult on the proton charge radius. This result is in a strong 
contradiction with a recent electron-proton scattering re- 
sult from MAMI [2 and the CODATA-2006 value Q, 
which basically originates from the hydrogen and deu- 
terium spectroscopy and involves large amount of exper- 
imental data and theoretical calculations. 

The discrepancy is at the level of 0.3 meV in terms of 
the muonic-hydrogen Lamb shift. Meanwhile the theo- 
retical uncertainty is equal to 0.004 meV and that from 
experiment is 0.003 meV. It is highly unlikely that the 
problem lies in either theory or experiment on muonic 
hydrogen. Nevertheless, it is important to clarify the 
theory of the muonic-hydrogen Lamb shift. 

We expect that the controversy will be resolved and 
the muonic-hydrogen Lamb shift will become the most 
accurate way to determine a value of the proton charge 
radius. For this reason it is important to have a reliable 
theoretical expression at the level of 0.003 meV [H. 

The theoretical expression consists of quantum- 
electrodynamics contributions and finite-nuclear-size cor- 
rections. 

After a calculation of all the light-by-light contribu- 
tions in order a^m^ 0, the quantum electrodynamics 
theory at this level of uncertainty is complete in a sense 
that all corrections have been calculated by at least one 
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author or one group. 

However, verification is required and we consider cer- 
tain corrections as not well established. That in partic- 
ular includes a contribution of the recoil effects in order 
a{Za)*mfj_ due to electronic-vacuum-polarization effects. 

The electronic- vacuum-polarization (eVP) effects form 
one of the most important features of a muonic atom 
which distinguish it from an ordinary atom. The lead- 
ing eVP efi^ect is due to a so-called Uehling potential 
and it produces the leading contribution to the Lamb 
shift in light muonic atoms. The correction is of order 
a(Za)^m^. Since that is the largest contribution, it is 
important to calculate the eVP terms including various 
higher-order effects. 

Another important feature is that the ratio of the mass 
m of the orbiting particle, a muon, is smaller than the 
nuclear mass M , but not so small as in an conventional 
atom and, in particular, m/M ~ 1/2000 in ordinary hy- 
drogen, while m/M ~ 1/9 in muonic hydrogen. That is 
why we need to find recoil corrections for most contribu- 
tions of interest. 

The purpose of this work is to obtain an a(Za)'^m^ 
contribution in all orders of m/M in light muonic atoms 
[Z = 1,2). The most straightforward way is to derive 
a non-relativistic expansion for both particles within a 
Breit-type equation 0| . However, for comparison we 
also use a Grotch-type technique [1] . The latter produces 
a non-recoil term and the m/M leading recoil correction, 
while the [m/M^ term is to be calculated separately. 
Both methods are considered in detail in this paper. 

These two methods are quite different. The Breit-type 
evaluation involves a non-relativistic expansion and the 
related perturbation theory includes terms of first and 
second order. Most of the first-order terms appear natu- 
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rally in momentum space, which is more challenging for 
a high-accuracy calculation. 

The Grotch-type approach involves non-trivial analytic 
transformations and as a result the expression for most of 
the energy contributions is obtained analytically in terms 
of a solution of the ordinary Dirac equation with a static 
potential. The only additional correction is expressed as 
a first-order perturbation, which for our purposes may 
be calculated using non-relativistic Coulomb wave func- 
tions. 

The results originally published for the relativistic re- 
coil contribution within these different approaches [9l-[Tl| 
are not consistent. The size of discrepancy is comparable 
to the experimental uncertainty 

Here, we find that in fact calculations in different ap- 
proaches are inconsistent because different gauges were 
used and in one of the calculations certain additional 
terms should appear. Those terms originate from re- 
tardation effects and from essential two-photon contri- 
butions and after proper corrections we find that the two 
approaches are consistent and produce the same result 
for the a{Za)'^m{iTi/M) terms. 

Comparing different approaches and comparing our re- 
sults with the results of other authors we basically focus 
our attention on muonic hydrogen. In addition, in sum- 
mary sections we also discuss other Z = 1,2 two-body 
muonic atoms. Units in which h = c — 1 are adopted 
throughout the paper. 

The paper is organized as follows. We first consider the 
non-relativistic leading eVP contribution and the leading 
relativistic term in order a{Za)'^m. Then we discuss dif- 
ferent gauges to take into account eVP effects. The two 
gauges we choose are closely related to static potentials 
applied in (9l4T]|. Both gauges are defined as a certain 
modification of the Coulomb gauge. 

We apply the Breit-type and Grotch-type approaches 
to calculate relativistic recoil contributions in the gauges 
and obtain consistent results in order a{Za)'^m{m/M). 
We find additional contributions due to a one-photon re- 
tardation contribution and due to two-photon exchange 
effects in one of those gauges. We demonstrate that the 
effective potentials generated by those additional terms 
agree with the difference in static potentials in the two 
gauges. 

We consider the term of order a{Za)'^m{m/ M)'^ only 
in the Breit-type approach. Since the value of this con- 
tribution depends on the definition of the nuclear radius, 
which may 'absorb' part of the correction, we recalcu- 
late a nuclear-finite-size term and obtain a semi-analytic 
result for it. 



II. THE LEADING NON-RELATIVISTIC AND 
RELATIVISTIC eVP TERMS AND THE 
RELATIVISTIC RECOIL EFFECTS 



The non-relativistic (NR) Uehling 
ily calculated both analytically [H, [k 



term can be eas- 
I and numerically 



(see, e.g., [Tj|) for an arbitrary hydrogenic state (with the 
reduced- mass corrections included). In particular, the re- 
sults for the low states in light muonic atoms are [l3 
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mn is the reduced mass 



mM 
M + m 



m and M are masses of the muon and nucleus respec- 
tively. We remind that in muonic hydrogen k ~ 1.5 and a 
characteristic value for the 2s and 2p states is k/2 ~ 0.75. 
In muonic helium ion this value is k/2 ~ 1.5. 

It is convenient to present the relativistic eVP correc- 
tion in order a(Za)^m as an expansion in powers of m/M 
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M 
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In this notation i5yp and Cq{k) are related to the 
leading relativistic correction, which, e.g., may be ob- 
tained from the Dirac equation including the Coulomb 
and Uehling potentials for a muon with the reduced 
mass. Such relativistic non-recoil contributions can be 
found through the Dirac-Uehling e quat ion both semi- 
analytically [ll HI] and numerically |10l . [l3 . [iTj for var- 
ious states. 
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In particular, for the n = 1, 2 states in muonic hydro- 
gen (k — 1.356146 . . . ) wc find from [16|] 

Co(ls) = -0.24488... . 

Co (2s) = -0.042224... . 

Co(2pi/2) = -0.0089077... . 

Co(2p3/2) = -0.00089011 . . . , (5) 

which indeed agree with the results from [l^, [HI, [13, [13 ■ 

The next coefficient, Ci, can be obtained in quite dif- 
ferent approaches (cf. |9l-[ll|) and one has to be care- 
ful while classifying contributions by counting photon 
exchanges. A rigorous consideration could start with 
a Bethe-Salpeter equation and, through rearranging its 
kernel, arrive at an effective one-particle or two-particle 
equation. When we refer here to one-photon exchange, 
we mean that the kernel of such an equation contains 
only a one-photon part. 

Any solution of an unperturbed problem is a sum- 
mation over an infinite number of such one-photon ex- 
changes. But all those many-photon exchange diagrams 
are reducible. 

The approach based on a Grotch-type equation (see, 
e.g., [l3l) immediately produces an appropriate summa- 
tion for the pure Coulomb exchange for the {ZaYm and 
{Za)'^m'^ /M terms. Meanwhile to take into account the 
eVP contribution one must use a perturbation theory, 
but only in its leading order. 

When the same physical contributions are treated non- 
relativistically, e.g., by applying the Breit-type approach 
(see d El), one has to use a certain effective (w/c)^ 
expansion. A perturbation theory is to be introduced 
already for a pure Coulomb problem to reach (Za)^m 
(both for the leading non-recoil term and for recoil cor- 
rections). As a result, a part of the relativistic one- 
photon contributions for a(Za)^m turns into reducible 
two-photon exchange contributions. 

One should distinguish between such reducible two- 
photon contributions and irreducible two-photon dia- 
grams which appear while calculating {Za)'^iv? /M in, 
e.g., covariant gauges. It happens that neglecting such ir- 
reducible contributions in one of the former calculations, 
namely in [To| , produces a result which is incomplete and 
must be corrected as considered below. 

The leading relativistic-recoil correction beyond the 
Dirac-Uehling term is of order a{Za)'^'W?' /M and may 
be calculated, as we mentioned, via various effective two- 
body techniques. 

The a{Za)^m term including various recoil correc- 
tions was calculated previously by a number of authors. 
Their results for E^p\2p - 2s) are 0.0169 meV [Ij, 
0.0169 meV ^ and 0.018759 meV [ll] and look con- 
sistent at first glance. 

However, as we note, the Dirac term 

- 2si/2) ^ 0.020843 . . . meV (6) 

has never been a problem and had been known from cal- 
culation of various authors for a while before the results 



mentioned above were achieved (see, e.g. [i3)- Once 
we subtract the Dirac term the results become strongly 
contradicting and the discrepancy is compatible with the 
uncertainty of 0.003 meV [J. 

In particular, for a value of -Eyp' — Eyp in the 
2pi/2 — 2si/2 splitting, the original results of the papers 
mentioned read -0.0041 meV in -0.0041 meV in ^ 
and -0.002084 meV in [III- 

The results obviously strongly disagree. We note 
that in principle the quoted calculations treated the 
{Za)'^m{m/M)'^ term differently. However, the later is 
smaller by a factor of m/M ^0.1 and cannot be re- 
sponsible for such a large difference by factor of two for 

p(rol) _ p(0) 
-^VP ^VP ■ 

Once we consider two-body diagrams, we have to real- 
ize that in principle there may be one-photon and two- 
photon contributions and in principle the one-photon 
term (see. Fig.[T]) includes a static part (found at fco — 0) 
and a retardation part (proportional to fco or fcp). 



X ^ X X ^ X 

(a) (b) 

FIG. 1: One-photon exchange diagram: for the pure photon 
exchange (a) and for the eVP contributions (b). 

The partial results, such as the non-recoil one-photon 
contribution, are not gauge invariant and only the sum 
of static one-photon, retardation one-photon and two- 
photon contributions is gauge invariant. 

The dominant part in any reasonable gauge is due to 
static one-photon exchange in order a{Za)'^m. (In prin- 
ciple, one can choose a ridiculous gauge with, e.g., a longi- 
tudinal part with a parameter 3> 1. That would allow to 
obtain a 'large' contribution beyond the mentioned above 
terms, but that does not have much physical sense.) The 
static one-photon exchange easily allows an efficient eval- 
uation both within the Breit-type and Grotch-type ap- 
proaches. 

However, those terms have already been covered by 
a consideration of the Dirac equation with the reduced 
mass. The purpose of this paper (as well as of (9l-[TT|) is 
to calculate recoil corrections in order a{ZaY'w?' /M and 
a{ZaY'm? /M"^ and we need to go beyond the leading 
terms. 

There is no 'objective' separation between static one- 
photon, retardation one-photon and two-photon contri- 
butions. E.g., considering the (Za)'^m contribution (in 
all orders in m/M) in different gauges, we have a kind of 
interplay of such terms. 
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Using the Coulomb gauge we find that for {Za)'^m 
there is neither a retardation correction nor a two-photon 
correction in order of interest and thus the complete re- 
sult in the {Zoif'm contribution may be achieved through 
an application of either the Breit equation (in all orders 
in m/M) or the Grotch equation for the leading m/M 
correction by applying the static one-photon kernel. 

Now we return to the results from |9l-[ll| cited above. 
All of them are obtained by calculating the static one- 
photon contributions. However, we demonstrate below 
that different gauges were in fact applied and only in one 
of them the retardation and two-photon contributions 
vanish while in the other they do not. 
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FIG. 2: Two-photon exchange diagram: for the pure photon 
exchange (a) and for the eVP contributions (b). 

We also note that the physical derivation of the Breit- 
type and Grotch-type equations should actually start 
from a two-body Bethe-Salpeter equation, next the lat- 
ter should be reduced to an effective one-body Dirac or 
two-body Schrodinger equation and afterwards contribu- 
tions of retardation and two-photon effects should be es- 
timated. The crucial part is to reduce the contribution of 
interest to an one-photon contjribution to the kernel of the 
effective equation (see, e.g., [18]). A further mathemati- 
cal transformation of the static one-photon contribution 
is rather a technical issue. 

A naive two-photon contribution, free of eVP, (see 
Fig.[2la) has infrared divergencies/singularities which in- 
dicate that it includes in fact a correction of lower order, 
namely the Coulomb correction (Za)^m. (A divergence 
appears once we neglect the atomic energy and momen- 
tum and consider the related diagrams as free scattering 
diagrams, otherwise we should speak about singularities.) 

The derivation through various effective approaches 
generates subtracted two-photon graphs (see, e.g., [11] )■ 
That is not a trivial issue. 

For example, if we choose the external-field approach 
for the first approximation, then we should some- 
how 'upgrade' {Zaf'm up to {Zaf'mji. The missing 
{Zafm^/M term comes from a two-photon contribu- 
tion. This example shows that a rearrangement of the 
diagrams applying a proper subtraction is crucially im- 
portant. 



The {ZaYm^ /M correction is a result of the calcula- 
tion of the nuclear-pole contribution of the two-photon 
box diagram. The effective-Dirac-equation approach sug- 
gests a complete subtraction of the pole of the heavy par- 
ticle JJJ (see Fig.[3]for the pole structure of unsubtracted 
two-photon diagrams of Fig[2]3; the pole of interest is de- 
noted as N-). 
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FIG. 3: Poles for the unsubtracted two-photon exchange dia- 
grams (see Fig. [2)). The upper plot is for the ladder diagrams 
(right graphs in Fig. [Sjl , while the lower one is for the cross 
diagrams (left graphs). For muonic hydrogen with k ~ 1 the 
pole structure for free diagrams and for eVP contributions 
(see Fig. [2] a and b, respectively) is similar. Here, A'' stands 
for nuclear poles, 7 is for photonic poles and /i is for muonic 
poles. 

Once such a pole is subtracted, we sec that the one- 
photon contribution is the only contribution for [Za)'^m 
terms in all orders of m/M as well as the non-recoil part 
of the (Za)'*m contribution. However, a recoil part of the 
{Za)'^m term arises in different gauges in different ways. 
The situation for the eVP (see Fig O*) is quite similar 
and the pole structure is also similar. 

Because of this similarity we briefly recall the situation 
with the {Za)'^m terms (in all orders in m/M) in different 
gauges. 

The electromagnetic interaction of two particles in dif- 
ferent gauges is determined by the shape of the photon 
propagator D^^{k). The term {Za)'^m originates from 
the Z?oo(fc) component, while {Za)^m corrections come 
from all components of Df^,j{k). 

One can immediately see that in any covariant gauge 
(in contrast to the Coulomb gauge) the Dqo compo- 
nent depends on the energy transfer fco and the retar- 
dation one-photon exchange produces a correction of or- 
der (Za;)^TO(m/M)^. In the case of non-zero values for 
the DiQ components certain terms of order {Za)'^{m^ /M) 
can also appear from related one-photon contributions. 

The static one-photon contribution obviously differs in 
different gauges and, after taking into account the retar- 
dation terms, the one-photon contributions still differ. 
One should take into account the two-photon diagrams 
to obtain complete {Z a)"^ (m^ / M) and (Za)'*m(m/M)^ 
contributions, which are indeed gauge invariant. 

It is easy to estimate a 'nominal' order of a two-photon 
diagram suggesting that it converges if we neglect all the 
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atomic effects. The order is (Zctf'm} jM. To obtain a 
lower order in Za, such as (Za)^m^/Af, we have to find 
terms divergent at low momentum. 

After the heavy-pole contribution is subtracted com- 
pletely, the only potentially divergent contributions are 
due to photon poles, if we close the contour in the lower 
half-plane (see Fig. [3]). 

The two-photon contribution in the Coulomb gauge 
has only a logarithmic divergence, which cannot change 
the fact that two-photon effects contribute in order 

If we consider another gauge, such as the Feynman or 
Landau gauge, the (Za)^m^/M terms do appear from 
the photonic pole contributions. We can see that such 
poles are important only in diagrams with Z?oo compo- 
nents for both photons, or with Dio{q) contributions. 

In the Coulomb gauge Dio = and the Dqq component 
of the photon propagator does not produce a pole and 
technically that is why the two-photon exchange in the 
Coulomb gauge does not produce any (Za)^(m^/Af) and 
(Za)*m(m/M)^ contributions. 

Once the heavy-particle pole is subtracted, the one- 
photon contribution is the only contribution for the 
{Za)'^m terms in all orders of m/M in the Coulomb 
gauge. 



III. eVP-CORRECTED PHOTON 
PROPAGATOR IN VARIOUS GAUGES 

Taking into account the transverse structure of the 
photon self-energy caused by the vacuum polarization 
tensor 



(7) 



one can derive in the Landau gauge 



Indeed, neither the Landau nor Feynman gauge is well- 
suited for the bound-state calculations and it is helpful 
first to perform a gauge transformation on the photon 
propagator ^ to reach a more suitable gauge. 

The photon propagator with an eVP correction in an 
arbitrary gauge can be presented in the form 



(13) 



where x — x{k) is an arbitrary function of k. 

We expect that in recent calculations [9l-[TTj| of the rel- 
ativistic recoil corrections to eVP different gauges were 
used. In this paper we consider two choices of the gauge 
function xi^) and two gauges. 

The gauge function can be presented as an expansion 
in powers of a 



TT 



(1) 
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For the x'''^' we chose the transformation which is to pro- 
duce the Coulomb gauge (for the free propagator), while 
for x^^-* we consider two options, which are presented be- 
low. 

The first choice is 



dvpe 







2(fc2-A2)(k2 + A2) ' 
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P'' ('^) ^ p 
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V{k 



2M 
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fc2 



V{k^) 



(8) 



(9) 



which differs only by longitudinal terms resulting from 
an obvious gauge transformation. 

Here, for the eVP we apply Schwinger's parametriza- 
tion 



P(fc2)-- / dvp,{v) 



' A:2 - A2 



where 



A^ 



l>2(l-i;2/3) 

l-i;2 
4m2 



1-^2 



(10) 

(11) 
(12) 



Let us refer to it as the CleVP gauge. 
A second possibility we consider here is 

(1) ^ ^0 

= - j, ^"^- 2(fc2-A2)k^ ' 



dvpe 



hi 



2(A:2 _ A2)k2 ' 
a f\ 1 P 



dvpc 



1 



k^ - A2 



X <5. 



/c^ hj 
~k2" 



(16) 



We refer to it as the C2eVP gauge. 

In the static regime (i.e., fcg = 0) the choice in Eqs. 
reproduces the potential applied in [ll[, while the 
choice in Eqs. leads to a potential considered in [lo|. 
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We note that the gauge (fT5)) is similar to the Coulomb 
gauge in a sense that the Dqq has no dependence on kg 
and there is no DiQ component. That means that the 
static one-photon contribution should produce a com- 
plete result. There are two contradictory results for the 
2p~2s Lamb splitting in literature [9|,[ll| and we confirm 
the latter result. Our result for muonic hydrogen is^ 



1/2) 



O.Of 87589 meV. 



(17) 



Considering the gauge ([T6| we note that, while Dio — 
0, the Z^oo component of the propagator depends on fco 
through the eVP tensor in Eq. pO| and one not only 
has to check the static one-photon term, but also cal- 
culate the one-photon retardation part and examine the 
photonic pole contributions for the two-photon diagrams 
(see Fig.^). 

To check the consistency of the result obtained and 
to compare with other existing calculations, we perform 
below four separate calculations applying either a CleVP 
gauge or C2eVP gauge within either the Breit-type or 
Grotch-type approach. 



IV. CALCULATION IN THE CleVP 
GAUGE (HID 

Now, let us perform the calculations in the CleVP 
gauge ([T5t , which because of lack of retardation and two- 
photon contributions in order a{Za)'^m (in all orders in 
(Za)) should produce a correct result in an easier way. 



A. Breit-type calculation 



Here and in further considerations we ignore the nuclear 
spin terms, assuming that the results are for the center 
of gravity of the related hyperfine multiplet. 

The result is obtained in the Coulomb gauge which is 
consistent with both CleVP and C2eVP gauge we are 
to consider. Here the first term is responsible for the 
(Za)^TO terms and the second produces {Za)'^m contri- 
butions (in all orders in m/M). To find the former one 
has to solve the related Schrodinger equation (let us de- 
note the energy and wave functions as Eq and ^'o) and 
to find the latter one has to find a matrix element of Hi 
over the Schrodinger-equation wave functions 4*0. 

Applying the CleVP gauge we find the additional 
terms in momentum space which are necessary to take 
into account eVP effects 



^cVP ^ H'^'^ + Hf^'^ 



cVP 



HI 



cVP 



Vuik) = -Aa{Za) / dv p. 



-4a(Za) / dv p, 



1 




4Mm 



(Pf 



(k2 + A2) 



(19) 



We start with a Breit-type calculations. 

First, we note that as it is well known (see, e.g., Q) the 
energy of hydrogenic levels without eVP can be found by 
considering a non-relativistic Schrodinger-type equation 
with a Hamiltonian 

H = Hq + Hi , 



Ho 



2m I 



Vc = 



p2 Za 

2m /J r 



8 \m^ AP 2 V 



ZOLTX / I 1 \ „T , ^ 



ZOL 



(r • (r •p)p) 



2mM r \ r 
Za f 1 1 
V 2mM 



((T-[rxp]) 



(18) 



^ Unless otherwise stated, the uncertainty is equal to unity in the 
last presented digit. 



The related expressions in coordinate space are 



:VP 



= Vu{r) = (Za) / dv p, 



HI 



gVP 



1 



1 



8m2 8Af2 

V4m2 
1 . 



2mM 



2mM 
1 

2mM 



V r>2 

r I 



(20) 



The Hamiltonian iJ°VP completely agrees with the one 
appearing in Q. 

Again, the first term is responsible for the a[Za)^m 
terms and the second one produces a{Za)'^m contribu- 
tions (in all orders in m/M). However, now the proce- 
dure is somewhat different. We are interested only in the 
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first order in a results and thus we can consider both where 
Hq^^ and Hf^^ as perturbations. 



To find the leading non-relativistic terms (see ([T])), 
we have to calculate (5'o|i?o^^|^'o), however, a simi- 
lar contribution of higher order, (5'o|-ffi^^|^'o), gives 
only a part of the a(Za)'*rn result. The other part 
results from second-order perturbation theory on the 
Schrodinger equation and it is of the form (see, e.g., [1]) 



2(vI/„|iJ§VP 



1 



(E - Ho) 



7^1 l*o) 



where the reduced Green function 



(21) 



is applied. In 



(E-Ho) 

our calculations for the non-relativistic reduced Coulomb 
Green function we used its presentation in terms of 
smaller and larger radii (cf. |17t). 

The wave function Vfo is expressed in terms of the re- 
duced mass mji while the dependence on m and M is due 
to apparent factors in Eqs. (ITS)) and (|20l) which allows to 
express results of our calculation of various matrix ele- 
ments in terms of Co, Ci and C2 as defined in Eq. 
That is helpful for further comparison with other calcu- 
lations. 

Indeed, the Co results reproduce the values ([5]) , as they 
should, while for other coefficients we find 



and 



Ci(ls) 
C\{2s) 
Ci(2pi/2) 

Cl(2p3/2) 



C2M 
C2(2.s) 

C2(2pi/2) 

^2(2^3/2) 



= 0.18153... , 
= 0.038089... 
= 0.00090127. 
= 0.00090127. 



-0.3631... 
-0.07618.. 
0.003542 . . 
-0.004475 . 



(22) 



(23) 



As we mentioned, the result for Ai?oVp(2pi/2 ~ '2^1/2) is 
consistent with the result of 111. 



K = cx ■ p + Pm 
1 



2M 



V 



(25) 



and 



To find a solution it is helpful, following "s], to intro- 
duce an auxiliary Hamiltonian 



Ki^ a-p + I3vi + V 



1 - pm/M 
1 - (m/M)2 



we note that 



K ^Ko + AK + O 



M' 



where 



2M 



4M 



and 



AK 



1 



2M 4M 



(26) 



(27) 



(28) 



(29) 



As it is known Q for the case of pure Coulomb poten- 
tial, AK = 0. Let us for the moment neglect /S.K for an 
arbitrary potential and look for a wave function of the 
form 

V'o=iV (^l-^[p',W-]) (1 + /?m)^. (30) 

Since the Grotch-type approach does not control the 
(m/M)2 terms, below we expand in m/M and neglect 
higher-order terms everywhere where it is possible. The 
results of such an expansion are denoted with 

In particular, we find for the normalization constant 



B. Grotch-type calculation 

The Grotch type of approach includes a few operations 
(see [sl for more detail). First, we have to present the 
two-body wave function as a product of a free-spinor for 
the nucleus and a four-component muon wave function. 
The potential is averaged over the nuclear spinor. 

The approach allows to reproduce the Dirac equation 
(with the reduced mass) and obtain the leading recoil cor- 
rection in order m/M. After we neglect terms of higher 
order in m/M ^ we arrive at an equation 



(24) 



N = 



1 



1 + 2^E/m + ^2 



E 
2M 



where the involved parameters are defined below. 

The final Grotch-type equation takes the form of an 
effective Dirac equation 



(a ■ p + /3m + v'j ip , 



V 



V 



V. 



(31) 



(32) 
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Solutions of the Dirac equation ((3T|) . and E, can be 
found, since the final equation takes form of a Dirac equa- 
tion with potential V and various effective parameters. 
The identities for "0 and E are of the same functional 
form as for a solution of a Dirac equation, but they ex- 
press the wave functions and energy in terms of effective 
parameters defined as 



E 
fh 




2M 

El - m'^/M 

™(l-t) 



/7"2 



El ~ m^/M , 



m 1 



1 



El 

M 



(33) 



Af2 



One can solve the equation (|5T|) as far as the solution of 
the conventional Dirac equation is known for a potential 
V with an appropriate accuracy. 

The final energy has a correction due to AK, which 
was neglected in order to obtain a solvable equation. The 
final result for the energy is 



E ^ Eo + {iPo\^K\ibo) ^ 



(34) 



We remind that AK cx ra/M and here we neglect all 
{m/MY corrections. 

Now, one can introduce the potential. In the C2eVP 
gauge, because the photon propagator is proportional to 
the free propagator in the Coulomb gauge, the equation 
would take the same form as for the Coulomb potential, 
but now with potential. 



V ^Vc + Vu. 



(35) 



However, in the case of the CleVP gauge, we have to 
introduce a certain correction, namely by redefining W , 



V 
W 

Wc 
Wu 



Vc + Vu. 
Wc + Wu, 
2Vc{q) 



%a{Za) I dv 



Pe{v) 



(g2 + A2)2 



(36) 



In both cases, we have to solve an effective Dirac equa- 
tion (pSj) with a potential (I55t . The energy levels for such 
a Dirac equation (with a reduced mass for the particle) 
with V — Vc are well known 

£;coui(n/j) = m + 4^„«^+4°j„, 



2n2 



2n3 



3 

4n 



1 



J + 1/2 



(37) 



and the linear in a correction due to eVP was calculated 
for the Dirac wave functions as explained in Sect. HIl (see 
for details Eqs. (P), (g]) and ^) 



Eyp — E, 



_ p(NR) 



VP 



E, 



(0) 
VP • 



(38) 



Solving the above equations, one can arrive at (see [iJ 
for details) 



Eq ~ Eq 

oul 

{nlj) + Eyp 



E. 



(NR) 
Coul 

M 



E. 



(NR) 
VP 



OK, ^ 



(39) 



We can now return to the AK term. The related cor- 
rection in the first order in eVP is a matrix element of 



AK 



1 

2M 
1 



1 



'iVuVc + -[VcAv\Wu]] 



(40) 



After a simple estimation of the operator, we find that 
it is sufficient to calculate the matrix element using 
Schrodinger-Coulomb wave functions, which are indeed 
well known. 

More detail on the application of the Crotch-type ap- 
proach to the eVP contribution will be published else- 
where [T9|. 

Finally, we obtain the results 

Ci(ls) ^ 0.18153... , 

Ci(2s) = 0.038087... , 

Ci(2pi/2) = 0.00090127... , 

Ci(2p3/2) = 0.00090127..., (41) 

which completely agree with the results ((22)) of the Breit- 
type calculation. 



V. CALCULATION IN THE C2eVP GAUGE JTeD 



For the C2eVP gauge p6)) we note that we have to 
calculate a static one-photon exchange, retardation one- 
photon contribution and two-photon contribution. Here, 
we first calculate the static one-photon contribution ap- 
plying the Breit-type and Grotch-type techniques and af- 
terwards we find the retardation one-photon contribution 
and two-photon contribution as a perturbation. 



A. Static one-photon exchange 

1. Breit-type calculation 

The Breit-type Hamiltonian is somewhat different from 
Eq. ([19]) and the addition is 



a{Za) r >^\Pt-Pf 
—n — I dv p, 
Mm / ' 



k2(k2 + A2)2 



(42) 
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which shifts the results for the static contribution. 
We find for the static one-photon contribution 



and 



Ci(ls) 
Ci(2pi/2) 

Cl(2p3/2) 



C2(ls) 

C2(2s) 

^^2(2^3/2) 



0.3981... 
0.06357... 
0.002845 . . 
0.002845 . . 



-0.7961... 
-0.1271... 

-0.0003444 . 
-0.008362 . . 



(43) 



which indeed does not coincide with Eqs. (p2|) and 
because such a contribution is not gauge invariant. 



2. Grotch-type calculation 



(44) 



failed to reproduce exactly the numerical results llOj for 
the correction neither by expressions given in [lOj nor 
by those presented in [23|, where further details of the 
calculation were given. Apparently, expressions of (loj 
still contain misprints. 

What is more important, the results we obtained in 
this section also disagree with the result fld\. The de- 
parture grows systematically between muonic hydrogen, 
deuterium and helium. (The results for the latter are 
presented in Sect. I VII] of our paper.) 

Meanwhile, we have discovered that our expression for 
the energy correction agrees with one presented in [2^ 
(see Eq. (116); it is also reproduced in Appendix A 
of [3). Our numerical results can be reproduced if we 
modify the erroneous expression for the term {^Qa) 
presented in Appendix A of 10]. 

Thus, we conclude that the result [l3| for the rela- 
tivistic recoil correction in light muonic atoms is unfor- 
tunately both incomplete because of lack of two-photon 
contributions and incorrect (because even the partial cal- 
culation contains a numerical error). 



A similar correction should be introduced into the 
Grotch-type approach. As we already mentioned, since 
the C2eVP gauge is proportional to the Coulomb gauge, 
we can use the same kind of equation as for the Coulomb 
gauge with a potential 

V = Vc + Vu 

and the W function defined within the same functional 
relation as for the Coulomb potential, namely as 



W{q) 



2V{q) 



The effective Dirac equation, which does not involve W, 
is indeed the same as in CleVP gauge, while the AK 
correction proportional to W is different. 

Proceeding similarly to described in Sect. lIVBI we ar- 
rive at 



Ci(ls) 
Ci(2s) 
Ci(2pi/2) 

Cl(2p3/2) 



0.39818... 
0.063585.. 
0.0028496 . 
0.0028496 . 



(45) 



which is consistent with Eq. and somewhat disagrees 

with fT(|. 

In particular, our result for Eyp'' — Eyp for the 2p— 2s 
splitting for the static contribution is —0.0042785 meV 
which is to be compared with —0.0041 meV in flO*!. 

The relativistic recoil eVP correction in light muonic 
atoms was calculated by Borie in [l3] in a way somewhat 
different from, but consistent with our treatment here of 
the Grotch-type approach in the gauge C2eVP. 

In the recent paper [10|[some minor corrections to the 
earlier papers j20| and |21| were introduced. Still, we 



B. Retardation one-photon exchange 

The retardation one-photon contribution and an essen- 
tial two-photon contribution can be calculated directly 
as a perturbation, since they are already related to effec- 
tive potentials which are smaller by a factor [Za)'^ than 
the non-relativistic contributions. It is also important 
that the second order perturbation term (see (|2T|) ) is the 
same as in CleVP gauge. That is because the free term 
of eVP propagator is the same for both gauges (which 
determines Hi) and the static limit of the eVP term of 
the propagator also does not change (which determines 
jjcVPy xhat says that only first-order perturbation the- 
ory terms are essential, otherwise the second-order terms 
similar to (l2Tt would appear. 

One can immediately find a related effective addition 
to the Hamiltonian due to retardation 



a{Za) 



dv Pe 



k2(k2 + A2)2 



or in coordinate space 

a{Za) 



retard 



dvp{v) 



where 



47rM2 

[{p^g(r)}-2p2Q(r)p2 



QW = T^-(2 + '^Se)|T^ 



(46) 



(47) 



(48) 



The results of direct calculations of Ci and C2 in 
muonic hydrogen are compiled in Tables |I] and |lTl respec- 
tively. The retardation effects under consideration are of 
order a{Za)'^m{m/MY and contribute only to C2. 
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Contribution 


Ci(ls) 


Ci(2s) 


Cl(2pi/2) 


Cl(2p3/2) 




Static 




0.39818 


0.063585 


0.0028496 


0.0028496 


Retardation 


















Two-photon 




-0.21654 


-0.025483 


-0.0019484 


-0.0019484 


Total 


0.18164 


0.038102 


0.0009012 


0.0009012 




TABLE I: Ci coefficients in 


tne y^ZQ V r J 


;auge. 






Contribution 


C2(ls) 


C2(2s) 


C2(2pi/2) 


C2(2p3/2) 






Static 




-0.7961 


-0.1271 


-0.0003444 


-0.008362 






Retardation 


0.2165 


0.02548 


0.001948 


0.001948 






Two-photon 


0.2165 


0.02548 


0.001948 


0.001948 






Total 


-0.3631 


-0.07618 


0.003552 


-0.004465 





50l) terms: 



TABLE II: C2 coefficients in the C2eVP gauge. 



C. Two-photon exchange 

For the two-photon exchange we have performed a cal- 
culation of the photon pole contributions. One has to 
carefully consider splitting of the contributions. We are 
interested in those which are singular at low momentum. 
Actually the photon-pole contribution is divergent in a 
formal sense at high momentum, but such a divergence 
is in order a{Za)^m'^ /M and thus is of a higher order. 
One has to separate properly the low-momentum and 
high-momentum contributions, and after that only the 
low-momentum one is of interest. The results are sum- 
marized in Tables HI and HIl respectively. 

We note that the sum of the results in the C2eVP gauge 
produces a result consistent with those in the CleVP 
gauge. We may also find an effective potential induced by 
a low-momentum contribution of the two-photon kernel. 
It is of the form 



27rM J (27r)3 J ^'■(q2 + A2)2 


in in 



q2 (k — q)" 



(q2 - k • q) 



(49) 



For the diagonal matrix elements we can replace the ex- 
pression (|49l) by an effective potential of the form 



a{Za) J dv pe- 


1 1 



k2(k2 + A2)2 
1 



MrriR 



(51) 



VI. FINITE-NUCLEAR-SIZE CORRECTIONS 

As we mentioned in the introduction, the definition 
of the nuclear radius for different spin values can pro- 
duce additional corrections in order {Za)'^m{m/M)'^ and 
a{Za)'^m{m/M)'^ . For this reason we consider in this pa- 
per the finite-nuclear-size (FNS) corrections. While the 
relation between the eVP recoil contributions and such 
corrections is considered in the next section, here we re- 
visit the FNS contributions with inclusion of eVP. 

Such contributions are known for a while [13, [12, ; 
basically numerically [2l| . Here we present semi-analytic 
results. 

We treat the FNS effects non-relativistically. In this 
case the FNS correction is of the form 



AE^ 



FNS 



6 

^■^FNS + ^-^FNS + ^-^FNS ' 



(52) 



where the wave function is the result of the Schrodinger 
equation with the reduced mass and potential Vc + Vjj 
and Rn is the rms nuclear charge radius. The result is 
not vanishing only for the s states. 
The leading term for the splitting 



a4ns(25 



(Za)^ 
12 



{R%ml()mR 



(53) 



is applied for a determination of the nuclear size radius 
from the experimental Lamb shift value and it is impor- 
tant to calculate corrections to it. 

The two corrections can be obtained in a way similar 
to a non-relativistic calculation of the leading eVP cor- 
rection to the hyperfine structure (cf. [TtI. [24 [25l| ) . The 
first term is 



A<s = ^^^(l*(0)P-|*c(0)P) 



(54) 



MniR J '^'"P'^ k2(k2 + A2)2 



(50) 



We see that at the end of the day the final effective 
Hamiltonians in both gauges are the same. The addition 
in the C2eVP gauge (|42p for the static term is eventu- 
ally cancelled out by the retardation and two-photon 



It is convenient to present the eVP correction in the 
form 



.^(1) _.^(o) |^(0)p-|vl^c(0)p 



(55) 



where we express the correction in terms of the leading 
term and the eVP correction to the wave function at the 
origin. The latter was studied in [25l - l28| . 
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Atom 






A 


AEfns 


fiR 


5.1975 


0.0170 


0.0110 


5.2254 


/iD 


6.0732 


0.0205 


0.0132 


6.1069 




102.52 


0.520 


0.323 


103.37 




105.32 


0.536 


0.333 


106.19 



TABLE III: The non-relativistic finite-nuclear-size corrections 
for tiie 2s state in light muonic atoms (in meV- (-Riv/fm)^). 
Here, fiHe stands for muonic helium ions. 



Atom 


Ai5oVp(2pi/2 — 2si/2) 


AScVp(2p3/2 - 2pi/2) 


AtH 


0.018759 


0.0049638 


/xD 


0.021781 


0.0057361 


M^He 


0.50934 


0.26920 


/i-'He 


0.52110 


0.27502 



TABLE IV: Relativistic recoil corrections to the Lamb shift 
and fine structure in light muonic atoms (in meV). 



The other term is 



6 



It was also found in \2t 



FNS 



(2s) = -G2.(/c/2)a4ns(2s) 



(56) 



(57) 



where 



G2s{z) - 



24 - 44z^ - 292-^ 



TT 

3z3 ' 36z2(l-z2)2 
- 20z2 + 33z4 - 20z^ + 8z 

2\2 



22z^ 



12z3(l-z2) 



■A{z), 



where A{z) is defined by Eq. ([2]). 

The numerical results on the finite-nuclear-size correc- 
tions for the 2s state in light muonic atoms are sum- 
marized in Table IIIII Note that for the Lamb shift the 
signs of the corrections are opposite. The results are 
slightly different but rather consistent with ones pre- 
sented in |lfl].2 



VII. RESULTS FOR THE LAMB SHIFT IN 
LIGHT MUONIC ATOMS 

To conclude a calculation of the relativistic recoil cor- 
rections, we have to fix the definition of the nuclear 
charge radius, which varies in literature. That is im- 
portant for the relativistic recoil corrections because with 
different definitions of the nuclear charge radius a certain 
part of the {Z a)'^ {m / M)'^m and a{Za)^{m/M)'^m terms 
may be incorporated in the nuclear-finite-size term. 

The nuclear spin takes different values for light muonic 
atoms, namely, / = 1/2 for a muonic hydrogen and 
muonic helium-3 ion, 7 = for a muonic helium-4 ion 
and 7=1 for muonic deuterium. The different nu- 
clear spin values are related to different effective two- 
body Breit-type equations for structureless particles (see, 
e.g., MM)- 



^ We are grateful to E. Borie who drew our attention to inaccura- 
cies in results for muonic helium. 



To be consistent with the experimental determination, 
we use the same definitions of the nuclear charge radius 
as applied in [l|, i, IMS ■ 

In this convention the Zitterbewegung term is present 
for half-integer spin nuclei, and not present for the integer 
case. 

The related calculation produces the results summa- 
rized in Table HVl 

Our results for the Lamb shift in light muonic atoms 
agree with the results of ill] and disagree with the results 
of [13 and m. 

For the fine structure, the results are also presented in 
Table lYl 

The result has been obtained in the CleVP gauge 
within Breit-type calculations and to control them we 
also performed a Grotch-type calculation in the same 
gauge. They are in perfect agreement with each other. 



VIII. SUMMARY 

Concluding, relativistic recoil contributions in orders 
a{Za)'^m{in/M) and a{Za)'^m{m/M)'^ to the Lamb 
shift in muonic hydrogen were revisited. The results 
published previously by various authors and obtained by 
different methods are inconsistent. In particular, result 
of the Breit-type calculation in [T^ is twice smaller than 
the related result from the Grotch-type evaluation in fl(| . 
The value of discrepancy, 0.002 meV, is comparable with 
the experimental uncertainty of 0.003 meV [l|. 

We perform here an evaluation in both approaches and 
find that the discrepancy between [ll| and [l3| is caused 
by the fact that both calculated the same value, namely 
the static one-photon exchange, which is not gauge in- 
variant. They applied different gauges. The gauge in- 
variant value is a sum of the static term and two other 
contributions, which are the retardation correction and 
two-photon contribution. While they are absent for the 
calculation of the a{Za)'^m{m/M) term in one gauge, 
they are not vanishing in the other. 

Once such contributions are taken into account and 
a relatively small numerical error in calculation [10] is 
fixed, we find perfect agreement between two approaches. 
Our results agree with those of |11|. 

We also consider other light atoms and perform cal- 
culations for the Lamb shift muonic deuterium and two 
isotopes of muonic helium. For a muonic helium-4 ion we 
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agree with [TJ] and disagree with [s^l- 

For control purposes we have also performed calcula- 
tions of the fine structure in order a{Za)'^m{m/M) and 
a{Za)'^m{m/M)'^. The resuh for the latter can be com- 
pletely restored from the result from Dirac equation (see, 
e.g., [16|). Our result for muonic hydrogen is in agree- 
ment with [s^j. 
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